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ABSTRACT 


The  theoretical  accuracy  and  coverage  of  the  two  most  promising 
air  traffic  surveillance  satellite  techniques  - -a  two-satellite  multi-rho 
technique  and  a  single  satellite  rho- theta -theta  technique- -are  determined 
by  the  application  of  a  mathematical  model  used  to  derive  error  isograms 
for  several  combinations  of  satellite  positions  and  range  and  angle  measure¬ 
ment  error  assumptions.  The  complete  derivation  and  discussion  of  the 
model  employed  is  included,  since  the  model  is  sufficiently  general  to 
apply  to  a  great  number  of  situations  involving  position  determination 
other  than  air  traffic  surveillance.  Some  conclusions  as  to  the  relative 
merits  of  the  two  techniques  studied  are  given,  and  recommendations  are 
made  for  future  activity. 


vii 


INTRODUCTION 


The  Systems  Research  and  Development  Service  has  the  responsi¬ 
bility  of  studying  innovations  in  operational  techniques  and  introducing 
new  technical  advances  for  future  implementation  in  the  air  traffic 
control  system  of  the  Federal  Aviation  Agency  and  is  making  studies 
of  Navigation,  Communication,  and  Surveillance  systems  for  use  on 
tht  North  Atlantic  Air  Route  for  the  post-J970  era.  A  report  by  the 
Communication  Group  of  the  University  of  Michigan  Institute  of  Science 
and  Technology  C  lJ  *  recommended  a  synchronous  satellite  communi¬ 
cation  system  for  Air  Traffic  Control  over  the  North  Atlantic  Air 
Route.  Of  all  the  systems  considered,  the  satellite  system  was  deemed 
most  reliable  and  economically  feasible  for  implementation  in  the 
post- 1070  time  period.  On  the  basis  of  this  recommendation.  Systems 
Research  and  Development  Service  has  been  making  studies  (in-the- 
house  and  under  contract)  pertinent  to  the  design  and  synthesis  of  a 
synchronous  satellite  communication  subsystem  for  over-ocean 
communications.  Preliminary  studies  of  propagation  and  coverage 
factors  show  that  the  VHF  (Very  High  Frequency)  Aeronautical  Communi¬ 
cation  Band  can  be  ueed,  and  that  a  tv'o-sateilite  synchronous  orbit 
configuration  is  preferable  to  three-,  four-  and  five -satellite  configura¬ 
tions  [2  J . 


In  a  concurrent  study  of  position  surveillance  systems  for  Air 
Traffic  Control  over  the  North  Atlantic  Ocean  by  the  Navigation  and 
Control  Laboratory,  Institute  of  Science  and  Technology  at  the 
University  of  Michigan,  the  recommendation  was  made  that  a  position 
data  acquisition  system  in  the  post-1970  time  period  consist  of  a 
two- satellite  synchronous  orbit  configuration  using  a  multiple  ranging 
technique.  Two  systems  embodying  this  technique  have  been  studied. 
The  first,  identified  by  the  acronym  LOCATES  (Location  of  Air 
Traffic  Enroute  by  Satellite;  was  proposed  by  the  University  of  Michi¬ 
gan  as  part  of  their  study  C31  .  The  second  was  a  two-satellite 
synchronous  altitude  systenr  for  Coverage  of  the  North  Atlantic  area, 
proposed  by  the  General  Electric  Company  as  an  alternate  to  the 
multi -satellite,  medium  altitude,  worldwide  system  of  their  final  report 
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All  references  will  be  listed  at  the  end  of  this  report. 


on  a  study  of  satellites  for  navigation  for  the  National  Aeronautics  and 
Space  Administration  C4^. 


The  Westinghouse  Electric  Corporation's  Defense  and  Space 
Center  under  Contract  NASw-785  sponsored  by  the  National  Aero¬ 
nautics  and  Space  Administration  has  proposed  a  navigation  satellite 
system  (which  can  be  used  for  a  position  data  acquisition  function) 
consisting  of  a  single  synchronous  satellite  using  a  combination  of 
ranging  and  angle  measurements  stationed  above  the  equation  between 
30°  and  40°  West  Longitude  C5]  . 

Although  eac*-  ~>f  tne  study  efforts  just  mentioned  presented 
some  data  on  the  positional  errors  to  be  expected  in  operational 
systems,  the  assumptions  made,  the  mathematical  methods  emp’oyed 
and  the  means  of  presenting  results  differed  for  each  study.  As  a 
consequence,  it  has  been  difficult  to  make  valid  comparisons  among 
the  systems.  Therefore,  an  in-house  task  was  established  to  construct 
a  general  mathematical  model  which  would  relate  the  distribution  of 
error  in  determining  the  position  of  aircraft  by  various  position 
finding  systems  to  the  distributions  of  assumed  input  errors  (bias 
errors  are  assumed  removable)  involved  in  measuring  angles  and 
ranges  between  station  and  aircraft  thereby  determining  the  theoretical 
system  quality  of  any  navigation  and  surveillance  systems  which  might 
be  hypothesized. 

The  body  of  this  report  consists  of  two  parts.  The  first  part 
(General  Theory)  includes  a  brief,  non-mathematical  discussion  of 
generic  concepts  ,  which  unify  position  determination  systems,  in 
terms  of  the  geometry  of  quadric  surfaces.  A  discussion  of  the 
problem  of  system  ambiguities  and  the  reasons  for  the  existence  of 
multiple  values  is  given.  A  general  Probabilistic  Position  Determina¬ 
tion  System  Model  is  then  constructed  and  developed  which  can  be 
applied  to  specific  navigation  and  surveillance  systems.  The  model 
is  a  first  approximation  to  the  real  physical  systems  under  considera¬ 
tion  whic'i  allows  using  the  tools  of  probability  and  statistics  in  their 
most  fundamental  aspects,  rather  than  the  more  esoteric  nuances  of 
these  mathematical  formalisms.  For  example,  it  is  assumed  that  the 
random  errors  are  normally  distributed  without  any  attempt  to 
justify  using  normal  statistics  by  introducing  the  Central  Limit  Theorem 
since  is  it  assumed  beyond  the  scope  of  the  report  to  do  so. 
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In  the  development  of  the  model,  coordinate  transformations  are 
made  by  introducing  matrix  methods  which  simplify  the  notation  and  allow 
a  clearer  picture  of  the  methodology  involved.  The  matrix  methods  also 
give  a  better  picture  of  the  transformation  of  axes  which  is  necessary  in 
a  model  to  be  used  for  computational  purposes. 

In  the  second  part  of  the  body  of  this  report  (Applications),  the 
mathematical  model  from  the  General  Theory  is  applied  first  to  the  two- 
satellite  system  (LOCATES)  and  then  to  the  one -satellite  system.  These 
two  systems  are  representative  of  air  traffic  surveillance  satellite  techniques 
being  studied  today.  This  section  includes  error  isograms  for  the  two 
systems,  compares  them  with  respect  to  positional  accuracy  and  coverage, 
and  makes  some  observations  of  the  relative  state-of-the-art  of  the  two 
sys  terns . 


GENERAL  THEORY 


1.  Fundamental  Concepts .  The  quantities  measured  by 
position  determination  systems  (directly  or  indirectly)  are  distances  and 
angles.  The  measurements  are  made  relative  to  known  or  measured 
positions  and  orientations.  The  generic  terms  RHO  and  THETA  are 
applied  to  distance  (range)  and  angle  measurements,  respectively,  and 
are  used  for  system  description.  For  example,  a  RHO-RHO-THETA 
system  measures  two  distances  and  one  angle.  From  each  measurement 
rr  ade  by  a  system,  a  quadric  surface  (a  surface  in  three-dimensional 
space)  is  generated.  A  position  in  three-dimensional  space  is  deter¬ 
mined  by  the  intersection  of  three  quadric  surfaces.  The  quadric  surfaces 
generated  by  measurements  of  position  determination  systems  are  limited 
to  spheres  generated  by  RHO  measurements;  cones  generated  by  THETA 
measurements;*  and  planes  generated  by  THETA  measurements.  (Strictly 
speaking,  the  plane  is  not  a  quadric  surface  since  the  equation  of  a  plane 
is  linear  but  the  linear  equation  is  a  degenerate  form  of  the  general 
equation  of  a  quadric  surface. )  When  a  position  is  to  be  determined, 
measurements  are  made  and  three  equations  are  given: 


M,  =  C, 
M2  =  C2 
M3  -  C3 


O) 


where  the  are  measures  of  either  ranges  or  angles  and  the  Cj  are 
constants.  The  equations  (1)  are  the  equations  of  three  quadric  surfaces 
in  their  most  general  form.  They  can  be  transformed  to  the  more  familiar 
Cartesian  coordinates  as: 

^The  cones  approximate  hyperboloids  for  some  systems. 
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F(x,  y,  z)  =  0 

G{x,  y,  z)  =  0 

H(x,  y,  z)  =  0 


U) 


where 


F (x,  y,  z)  =  a, ,  x2  +  2  a,  2  x  y  +  2  a, 3  x  z 

+  a22yJ  +  2  a23y  z  +  a33  z3 

+  2  al4  x  f  2  a24y  +  2  a34z  +  a44 
with  the  ajk  =  constants. 

G(x,  y,  z)  and  H(x,  y,  z)  have  the  same  form  but  the  constants  aik  are 
different.  In  order  to  determine  a  position  in  x,  y,  and  z  space.,  it  is 
necessary  to  solve  the  equations  (2)  simultaneously  for  x,,  y,,  z,,  the 
coordinates  of  the  position.  Because  the  equations  are  quadratic  (the  in¬ 
dependent  variables  are  second  degree)  solving  them  simultaneously  yields, 
in  general,  eight  positions.  For  specific  surfaces  and  a  particular  choice 
of  origin  (usually  the  center  of  system  symmetry)  they  are: 

P,  =  (x,.  y|t  zt) 

P2  =  (x,,  y,,  -z,) 

P3  =  (x,,  -y,,  z,) 

P4  =  (x„  -y,,  -z,) 

Ps  =  (-x„  y,,  z,) 

P*  =  (-x,,  y,,  -z,) 

P7  =  (-x|(  -y,,  z,) 

Ps  =  ("X|,  -y,,  —  z,) 


The  number  of  points  found  for  any  specific  position  determination  system 
depends  upon  the  surfaces  generated  by  the  system.  A  RHO-RHO-RHO 
system,  for  example,  generates  three  spheres.  The  intersection  of  two 
of  the  spheres  is  a  circle  and  the  intersection  of  the  circle  with  the  third 
sphere  yields  two  points.  Systems  can  be  engineered  to  determine  a 
unique  point  in  space.  For  example,  a  tracking  radar  (RHO- THETA- 
THETA)  measures  a  range  and  two  angles: 


R  =  R, 

0  =  0, 

<t>  =  $1 


(3) 


where  R  is  range,  0  is  elevation  angle,  $  is  azimuth  angle,  and  R,,  9,  , 

<X> , ,  are  constants.  The  equations  (3)  are  the  equations  of  a  sphere,  cone, 
and  plane,  respectively.  The  point  (R,  ,  9  ,  ,  <J>,  )  is  the  intersection  of 
the  quadric  surfaces  as  shown  in  Fig.  1. 
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FIG.  1  QUADRIC  SURFACES  ASSOCIATED  WITH  A 
TRACKING  RADAR  SHOWING  TH/SIR  INTER¬ 
SECTION  AT  (R|,  0,,  <$,). 


The  point  (R|,  0,,  <$,)  is  unique  because  of  system  conditions.  R  =  -Rt 
yields  the  same  sphere  as  R  =  R|  but  the  system  always  measures 
R  >  0.  Since  the  system  scans  the  half  space  above  earth,  it  only 
measures  0  0.  The  system  is  set  up  so  that: 

<t>  =  <X>i  *  <$,  +  160° 

i.e.,  it  has  a  unidirectional  beam  of  radiation.  Because  of  these 
auxiliary  conditions  the  tracking  radar  determines  a  unique  point.  If 
any  coordinate  transformations  are  used  on  the  tracking  radar  system 
equations,  the  auxiliary  conditions  must  be  carried  along  to  avoid  am¬ 
biguities.  The  problem  of  multiple  values  of  a  system  can  usually  be 
handled  quite  easily,  but  some  consideration  of  them  must  be  made  in 
any  system  design. 

2.  Probabilistic  System  Model.  It  has  been  shown  that  a  position 
in  space  is  determined  by  finding  the  intersection  of  three  quadric  sur¬ 
faces.  When  a  position  is  to  be  determined,  measurements  are  made 
giving  the  values  M,,  M2,  M}.  These  measurements  are  used  to 
generate  a  set  of  equations  which  can  be  solved  simultaneously  for  the 
position  coordinates.  M,,  M2,  and  M}  are  rarely  exact  measures  (the 
probability  that  they  are  exact  is  zero)  since  random  errors  can  never 
be  completely  excluded  from  them.  Since  we  are  interested  in  finding 
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the  probability  that  a  measured  position  is  within  some  linear  distance  of 
the  true  position,  we  assume  that  the  i'th  measurement  of  a  position  in¬ 
cludes  linear  errors  Cm,,  Cm^  and  CA3in  the  measures  of  M,,  M2,  and 
Mj,  respectively.  If  the  measurement  is  a  range  measurement,  the  error 
is  along  the  vector  from  the  station  to  the  position,  and  if  the  measure¬ 
ment  is  an  angle,  the  error  is  linear  and  normal  to  the  vector  from  the 
station  to  the  position.  The  angular  error  in  an  angle  measurement  can 
always  be  determined  by  dividing  the  linear  error  by  the  magnitude  of 
the  vector  from  the  station  to  the  position.  If  a  new  measurement  (the 
j'th)  is  made  on  the  system  with  no  change  in  the  actual  position,  a  new 
set  of  linear  errors  results  and C^3 where,  in  general. 


0*3*  01*3 


If  a  whole  series  of  measurements  is  m^de  with  the  actual  position  un¬ 
changed.  a  distribution  of  position  determinations  about  the  true  position 
is  given.  From  the  distribution  a  probability  density  function  is  de¬ 
rivable.  When  the  probability  density  function  for  a  region  of  space  i« 
known,  it  is  possible  to  determine  the  probability  that  a  position  deter  ¬ 
mination  will  fall  into  any  portion  of  the  region.  If,  for  example,  the 
probability  density  function  throughout  some  given  volume  of  space  V  is 
p,  the  probability  P  that  a  position  determination  is  in  V  is  given  by: 


P 


(4) 


The  model  to  be  constructed  in  this  report  can  be  applied  to  specific 
position  determination  systems  and  the  result  of  the  application  will  give 
the  probability  that  a  position  determination  is  in  some  given  region  of 
space  about  the  true  position. 

The  integration  of  equation  (4)  is  easiest  to  perform  when  the 
probability  density  function  p  and  the  volume  V  over  which  the  integration 
is  to  be  made  are  geometrically  compatible.  For  example,  if  p  is  a 
spherically  symmetrical  function,  and  the  integration  is  over  a  spherical 
volume  whose  center  is  the  center  of  symmetry,  the  problem  has  geo¬ 
metrical  compatibility  and  is  in  its  most  suitable  mathematical  form.  For 
a  system  where  errors  are  considered  random  and  normally  distri¬ 
buted,  the  natural  geometrical  configuration  is  the  ellipsoid.  That  is,  the 
surfaces  on  which  the  probability  density  function  is  constant  are  homo- 
thetic  ellipsoids  about  the  true  position  as  the  center  of  gravity  of  the 
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distribution  when  the  errors  have  zero  means.  Thus,  the  ellipsoids  are 
called  trivariate  normal  equi-probability  error  density  ellipsoids.  Al¬ 
though  ellipsoids  can  be  described  in  any  arbitrary  orthogonal  coordinate 
set,  the  equations  are  generally  complicated  and  difficult  to  work  with  so 
it  is  usually  best  to  work  in  the  coordinate  axes  which  are  coincident  with 
the  principal  axes  of  the  ellipsoid.  Since  many  real  position  determination 
systems  are  originally  given  in  the  natural  physical  set  of  coordinate  axes 
which  are  generally  non-orthogonal  and  not  coincident  with  the  principal 
axes  of  the  system  ellipsoids,  it  is  desirable  to  transform  from  the 
original  non-orthogonal  set  (M,,  M2,  M3)  to  the  orthogonal  principal  axes 
set  (p( ,  p2,  p3).  The  systems  whose  original  physical  sets  of  coordinate 
axes  do  coincide  with  the  principal  axes  sets  are  special  simple  cases 
included  in  the  general  mathematical  model  to  he  developed  here.  Since 
the  original  set  of  axes  is  assumed  non-orthogonal,  the  transformation 
to  principal  axes  is  made  in  two  steps.  First,  a  transformation  matrix, 
T.is  developed  which  transforms  errors  from  the  given  non-orthogonal 
set  (M,,  M2,  M3)  to  an  arbitrary,  but  known,  orthogonal  set  of  axes 
(x,,  x2,  Xj).  Second,  a  rotation  matrix ,  R,  is  developed  which  rotates 
errors  from  the  orthogonal  set  (x|f  x2,  x3)  to  errors  in  the  principal 
axes  set  (p , ,  p2,  p3).  In  mathematical  form: 


The  complete  transformation  is  given  by: 


cp=  RTCA 


Before  the  T  and  R  matrices  are  developed,  the  significance  of  having  the 
errors  in  the  principal  axes  coordinate  set  will  be  discussed.  As  has 
been  mentioned  previously,  ellipsoids  can  be  described  in  any  arbitrary 
orthogonal  coordinate  set  (x,,  x2,  x3).  However,  the  probability  density 
function  in  the  (x,,  x3,  x3)  set  is  given  by: 


P  (CX,-  CXj'  Cx^  = 


e“  2  Q  (  Cx,-Cx2  .£Xj  ) 
(2r)^5  ®x,  <Tx2ax3^|i?| 


where: 


Q<Cx,’CX2,  Cx3)  =  -j^Se'k‘5^k 
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2  J 

^2  1 

^22  ^2  3 

^3  1 

^3  2  ^3  3 

(i,  k  =  1,2,  3) 
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e..  --  - 


r#xk n  £  (cx,cxj 

'  k  1  =  1 


ik  '  “3C  3 n”  x.  *-  xL;  <j»K  -  »  2»  3\ 

U  *  k  / 


1 


(j  =  1.2,  3) 


e-  tc; 

It  is  evident  that  p(Cx,fIx2,£oc3)  is  a  complicated  expression  and  that  the 
integration  of  equation  (4), 


■Iff 


p  =  I  IJ  p(£x,f£xj,Cx3  )  d V, 
V 


would  be  quite  difficult  for  any  volume  V.  However,  if  the  error  com¬ 
ponents  are  rotated  to  the  principal  axes  set  (p)t  p2 ,  p3)  according  to: 


€'  =  R  C 


then  the  correlation  is: 


Then: 


‘jk 


0  k  *  j 

1  k  r  j 


(j,  k  =  1,2,3) 


(5) 


1*1* 


1  0  0 
0  1  0 
0  0  1 


=  1 


CP. 


Q<cp,.Cp1.Cp))  +  %- 

*1  *2 

1  ^  2 

(£'  ) 


(j  =  1,2,3) 


pj  “  rr  '^p, 

The  probability  density  function  in  the  principal  axes  set  becomes: 

.i&  ♦  ^  ♦  %) 

y^p  °p  °P  / 

p,Cp„CpI.Cp)),  a - J- 


Pi  p2  p> 
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a  relatively  simple  expression.  To  find  the  probability  that  a  position 

determination  is  in  an  ellipsoid  whose  axes  are  ca„  ,  cOL  ,  cOL,  where 

f  .  ,  Pi  ri  P3 

c  is  f  constant,  evaluate  the  integral: 


Let 


'  VPi  Pi  °i>J 


dz, 


dCp,  dCp.dC  p3 


”P. 

=  d_£E L 


°p3 


dz,  = 


dC 


El 


dz,  = 


dC 


El 


Then 


dz,  dz2  dz. 


Transforming  to  spherical  coordinates 


Then 


z,  =  r  sin  6  cos  <J> 
z2  =  r  sin  0  sin  <J> 
z3  =  r  cos  8 


Z(  +  A  +  A  =  r2 

dz,  dz2  dz3  =  r2  sin  0  drd  9  d<$ 


where 


r2  d  r 


r2  d  r 


9 


Integrate  by  parts  letting 


then 


u  =  r 


du  =  dr 


and  dv  =  re  2  dr 


and  v  =  -  e  2 


f  J  2 

/  |c 

vH 

udv  =  yr* 

V  1° 

/•c  r2 

2 

[  ~  T* 

y/TT 

1  —  c  e  + 

}  9 

Jo 

Pc  =  2 


(^/  2  dr  -  •*  *) 


where 


Pc  =  2  (Ac  -  cBc) 


c  r2 
'  T  dr 


vbrl  e 

H 

L J'.-i 

\nr) 


c  j7 
e  ^  dr 


vrrf 


o  r2 
e'7"  dr 


B. 


dr  -  0.5 


_ 

2 


Ac  and  Bc  are  tabulated  in  many  probability  and  statistics  books  [dj.  A 
very  short  table  of  Pc  values  is  given  in  Table  I. 


TABLE  I 

A  SHORT  TABLE  OF  VALUES  OF  PROBABILITY  OF  A  POSITION 
DETERMINATION  BEING  IN  A  ONE,  TWO,  OR  THREE  SIGMA 
ELLIPSOID  WHOSE  CENTER  OF  GRAVITY  IS  THE  TRUE  POSITION 
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c 

Ac 

c  Bc 

P 

c 

1 

.  341.34 

. 24197 

.  19874 

2 

. 47725 

.  10798 

. 73854 

3 

. 99865 

.01329 

.97072 

Table  I  shows  the  probability  that  a  position  determination  in  a  one, 
two,  or  three  sigma  ellipsoid  is  .  19874,  .73854,  and  .97072,  respectively. 
The  actual  physical  dimensions  of  the  principal  axes  of  the  ellipsoid  are 
given  by  co^,  c  cr^,  cc^  (with  c  =  1,  2,  3,  respectively)  where  the  values 
of  Op  are  to  be  determined  from: 

‘v-i-S;  ,cp/  = *• 2> 3) 

j  >«i  j 

and  Cl  (j  =  1,  2,  3)  are  determined  from: 

Pj 

Sj,  =  RT 

Since  the  orientations  of  the  (M,,  M2,  M3)  axes  are  given  and  the  T  and  R 
matrices  will  be  developed,  the  ellipsoid  orientation  will  also  be  deter¬ 
mined. 

In  order  to  find  the  T  matrix  in: 

r-  X  M 

it  is  necessary  to  find  the  vector  Ci  from  the  given  components  C**, 

^■*4*  Assume  that  the  errors  are  small  and  the  distance  between 

stations  at  which  measurements  are  made  and  the  region  where  position 
determination  are  to  be  made  are  large.  Under  this  assumption  spheres 
and  cones  can  be  considered  planes  to  a  first  approximation.  That  is, 
for  a  sphere  as  large  as  the  earth,  and  short  distances  on  its  surface, 
the  earth  can  be  considered  flat.  Since  the  approximation  can  yield  only 
conservative  values  the  only  difficulty  which  may  occur  is  the  introduction 
of  singularities  in  the  ellipsoid  dimensions.  That  is,  ellipsoids  of  in¬ 
finite  dimensions  may  result.  To  obviate  such  difficulties  it  is  necessary 
to  introduce  higher  order  approximations  for  specific  systems  at  points 
where  singularities  appear.  Such  higher  order  approximations  would  not 
do  for  the  entire  system  analysis  because  of  the  rule  that  as  a  mathematical 
model  is  made  more  realistic,  it  also  becomes  more  complicated  and 
difficult  to  handle.  For  the  problem  of  this  report,  each  possible  position 
determination  would  require  a  separate  calculation,  whereas  the  first  order 
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approximation  will  give  a  solution  in  closed  form  valid  everywhere  except 
at  points  where  singularities  appear.  At  such  points  higher  order  ap¬ 
proximations  can  be  made  for  meaningful  results. 


M3 


FIG.  2  THE  ERROR  VECTOR  €*  *  IS  GIVEN  BY  THE 

INTERSECTIONS  OF  PLANES  DRAWN  NORMAL 
TQ  M„  M2,  M3  FROM  ERROR  COMPONENTS 
^  /i2'  *3 

— ♦  . 

In  Fig.  2  the  error  vector  G^  is  found  by  drawing  planes  normal  to  the 
axes  at  G  fa  and  G  ^ 3.  ^  gives  the  measured  position  Pj 

relative  to  the  true  position  PQ.  It  should  be  noted  that  the  vector  com¬ 
position  of  Fig.  2  does  not  correspond  to  that  of  the  ordinary  vector 
composition  of  vector  geometry.  The  method  used  here  is  necessary  to 
correspond  to  physical  position  determination  systems.  The  axes  M,f 
M2,  M3  are  generally  non-orthogonal.  Assume  a  new  set  of  axes  (x',, 
x'r  x3  ),  where  x{  is  in  the  plane  of  M,  M3  and  normal  to  MJf  x'7  is  in 
the  plane  of  M2  M3  and  normal  to  M3,  and  x3  coincides  with  M3  as  shown 
in  Fig.  3. 


M, 


FIG.  3  A  SKETCH  SHOWING  THE  TRANSFORMATION  FROM 
THE  (M„  M2l  M3)  SET  OF  COORDINATES  TO  THE 
(x,'  ,  x\  ,  x!,  )  SET 
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The  angle  between  M,  and  Mj  is  <J> ,  the  angle  between  Mv  and  Mj  is  4', 
and  the  angle  between  xj  and  xj  is  0  as  shown  in  Fig.  3.  Mote  that 
the  (xj,  xj,  xj)  set  of  coordinates  is  also  non-orthogonal  since  0  is  not 
generally  90°. 

Fig.  4  shows  the  M,  M3  plane  and  the  planes  drawn  normal  to  M,and 
M3  at  G  mi  and  Cm3  intersecting  at  Pj.  A  plane  from  the  intersection  Pj 
normal  to  the  xj  axis  marks  off  C^i  on  that  axis. 


M3.  x  j 


FIG.  4  THE  PLANE  M,  M3IS  SHOWN  with  THE  PLANES  . 

DRAWN  NORMAL  TO  M,  AND  M3  AT  G£,  AND  G 
INTERSECTING  AT  p[ 

From  Fig.  4: 


ci;  =  G/i,  CSC  <t>  -  Ci3  cot  <t> 

=  G* 


Fig.  5  shows  the  M2  M3  plane  and  the  planes  drawn  normal  to  M2 
and  M3  at  Gm2  and  G^jir  tersecting  at  P  i.  A  plane  from  the  intersection 
of  Pi  normal  to  the  x2  axis  marks  off  G  icj  on  that  axis. 


FIG.  5  THE  PLANE  M2  M3  IS  SHOWN  WITH  THE  PLANES 
DRAWN  NORMAL  TO  M?  AND  M3  AT  C  ^  AND  G  * 
INTERSECTING  AT  P’j  2 
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From  Fig.  5: 

C* .  =  Cm2  CSC  f  -  d3  cot  i 

Glj  =  ^3 

Transforming  from  the  (M,,  M2,  M3)  coordinates  to  the  (xj,  x^  ,  x!,) 
coordinates  gives: 

Cl,  =  Cm,  csc<t>-  C^3  cot<J> 

xi 

G  =  C^csc  *  -  Gm3  cot  * 


To  transform  from  the  (xj,  x2,  x3)  coordinates  to  the  orthogonal  coordinates 
(x,,  x2,  x3)  consider  Fig,  6,  Fig.  6  shows  the  xj  x2  plane  and  the  planes 
drawn  normal  to  xj  and  x2  at  G^j  and  G^i  intersecting  at  P'V  .  A  plane 
drawn  through  P'V  and  normal  to  the  x2  axis  marks  offG]^  on  that  axis. 


FIG.  6  THE  PLANE  xj  x}  IS  SHOWN  WITH  THE  PLANES 
DRAWN  NORMAL  TO  xj  AND  x^  AT  Gj,i  AND C^i 
INTERSECTING  AT  P".'  '  2 

From  Fig.  6: 

C*  =  C1,  esc  0  -  £'  cote 

*2  Xj 


Therefore: 
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CX,  = 

Cm,  csc  $  " 

CM3COt 

CX2  = 

(C  M2  CSC  i  - 

C M 3  COt  \^)  CSC  0 

-  (Cm,  esc  <t> 

-  C m3 cot <J) )  COt0 

Cx3  = 

Gi3 

In  matrix  notation 


1 

X 

CSC  <J> 

0 

-  cot  <!> 

re*  1 

4, 

= 

-  CSC  <t>  cot  0 

CSC  ^  CSC  0 

COt(J>  cot  0  -  cot  ^  CSC 0 

f  i 
^-m2 

.4,. 

0 

0 

1 

The  T  matrix  is  therefore: 


T  = 


esc  $  0 

-  CSC  <J>  cot  0  C8C  41  CSC  0 


-  cot  <J> 

cot  <|>cot  0  -  cot  ^  CSC  0 
1 


and  it  transforms  errors  from  a  non-orthogonal  coordinate  set  (M,,  M2,  . 

M3)  to  an  arbitrary  (but  known)  orthogonal  coordinate  set  (xi,  x2,  x3). 

When  the  errors  are  found  in  an  orthogonal  set,  they  can  be  transfoimed 
to  any  other  orthogonal  set  by  using  a  rotation  matrix  R,  In  particular  it 
is  desirable  to  rotate  from  the  (xi,  x2,  x3)  set  to  the  principal  axes  set 
of  coordinates  (  p,,  p2,  p3).  The  rotation  matrix  R  is  not  unique.  The 
most  commonly  used  one  specifics  the  rotation  matrix  in  terms  of  the 
three  independent  parameters  called  the  Eulerian  Angles.  Since  the 
many  sources  in  the  literature  are  not  consistent,  the  R  matrix  will  be 
developed  in  detail  to  show  the  transformations  involved.  The  rotation 
R  is  broken  up  into  three  separate  rotations.  Beginning  with  the  (X|,  x2, 
x3)  set,  rotate  counterclockwise  through  an  angle  A  about  the  x3  axis,  as 
shown  in  Fig.  7. 


FIG. 


SHOWS  A  COUNTERCLOCKWISE  ROTATION  ABOUT 
THE  x,  AXIS  THROUGH  THE  ANGLE  A 
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In  matrix  notation 


The  rotation  is  given  by: 

“  m 

cos  A  sin  A  0 

A*  =  -sin  A  cos  A  0 

0  0  1 

Next  rotate  counterclockwise  about  the  xj  axis  the  angle  B  as  in  Fig.  8. 


FIG.  8  SHOWS  A  COUNTERCLOCKWISE  ROTATION  ABOUT 
THE  xj  AXIS  THROUGH  THE  ANGLE  B 
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The  rotation  is  given  by: 


B 


cos  B  sin  B 


0  -sin  B  cos  B, 


Finally  rotate  the  (x't',  x,',  x")  set  the  angle  C  counterclockwise  about 
the  X]  axis  as  in  Fig.  9. 


FIG.  9  SHOWS  A  COUNTERCLOCKWISE  ROTATION  ABOUT 
THE  x'3'  AXIS  THROUGH  THE  ANGLE  C 


i 

cos 

C 

sin  C 

o" 

1  ^ 

—  ~x 
u 

1 _ 

ck 

= 

-sin 

c 

cos  C 

0 

Ciu 

X2 

kJ 

.  0 

0 

1. 

r  * 

s-  - 

x3 

The  rotation  is  given 

by: 

- 

cos 

c 

sin  C 

0 

C'  = 

-sin 

c 

cos  C 

0 

0 

0 

1 
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Because  the  three-dimensional  rotation  matrices  form  a  non-abelian  group 
the  order  of  rotation  is  important  [7].  In  general: 

A' B'  *  B'A' 

The  distributive  law  holds,  however,  so  that: 

(A'B')(C')  =(A')(B,C')  =A'B'C’ 

The  rotation  matrix  P  is,  therefore: 

R  =  A'B'C' 


Multiplying  out  yields: 


cos 

C 

cos  A 

cos  C 

sin  A 

sin  C 

sin  B 

-cos 

B 

sin  A 

sin  C 

+  cos  B 

cos  A 

sin 

C 

-sin 

C 

cos  A 

—sin  C 

sin  A 

cos  C 

sin  B 

-cos 

B 

sin  A 

cos  C 

+  cos  B 

cos  A 

cos 

C 

sin 

B 

8  in  A 

-sin  B 

cos  A 

cos  B 

so  that 


fex  1 

p, 

-C- 

xi 

CP, 

= 

R 

ei 

2 

el 

cj. 

P'1 

L  *3J 

and 


The  components  of  error  in  the  principal  axes  coordinate  set,  therefore, 
can  be  determined  using  R  and  T.  The  variances  in  the  principal  axes 
set  are  defined  by: 


(j  =  1,2,3) 
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The  variances  are  found  as  functions  of  trM),  <XM2,  cr*^,  <£,  0,  A,  B,  C. 

<JM),  °M3  are  given.  <J> ,  \j/ ,  0  are  determined  by  the  geometry  of  the 

particular  physical  system  under  consideration.  A,  B,  C  are  determined 
from: 


_1 


<j,  k  =  1,  2, 3 j 

\j*i<  I 


The  <Tp.  are  defined  completely  and  the  physical  dimensions  of  the  error 
density^ ellipsoid  are  known.  The  spatial  orientation  of  the  ellipsoid  is 
also  known  in  terms  of  the  angles  A,  B,  and  C.  The  probability  that  a 
position  determination  falls  into  a  one,  two,  or  three  <f  ellipsoid  is 
.19874,  .73854,  an^..  97072,  respectively. 


There  is  an  alternative  method  for  solving  for  the  Op.  using  formal 
matrix  methods  given  in  some  texts  [s]  which  will  be  shosAi  here  in  order 
to  derive  a  useful  approximation.  Rather  than  finding  error  components 
in  the  (p  ,  p2,  p3)  axes,  the  error  components  are  determined  in  the 
orthogonal  set  (x,,  x2,  x3)  by  using  the  T  matrix.  Then  the  so-called 
covariance  matrix  Ax  is  developed  where: 


■ 

ax, 

*.2 

cr 

xi 

x2 

a, 

x3 

— 

e„ 

c r 

a 

a2 

CT 

*2 

X| 

*2 

x3 

$ 

O' 

<T 

P 

a 

crl 

x3 

X. 

32 

x3 

x2 

x3 

by  using  the  definitions: 


1 


N 


ji  k, 

j  *  k 


?jj  *  1  (3=  1.2.3) 

The  formal  matrix  methodology  for  diagonalizing  the  covariance  matrix 
is  the  similarity  transformation  [9]: 


Ap  =  R  Ax  Rt 

where  RT  is  the  transposed  R  matrix  and  where: 


"p, 

°,2 

ar_  o_ 

Pi  Pj 

^3°P,  ^3 

f>'°p, 

^Pl 

<T2 

P2 

^23  aPj  ^Pj 

0 

P 

a  a 

cr2 

.  'Si  P3 

Pi 

'32 

P3  P2 

P3 

n 


Since  the  off-diagonal  elements  of  a  diagonal  matrix  vanish;  that  is 


J.  k  -  1.2.3# 
j  *  *  ( 


the  angles  A,  B,  and  C  of  the  R  matrix  can  be  determined.  The  tf’p., 
therefore,  are  determined  and  the  ellipsoid  orientation  is  also  founA  from 
A,  B,  and  C.  It  is  evident  that  the  probability  of  finding  a  position  de¬ 
termination  in  a  sphere  containing  an  ellipsoid  is  greater  than  the  proba¬ 
bility  of  finding  it  in  the  ellipsoid.  For  example,  the  probability  of 
finding  a  position  determination  in  a  sphere  containing  a  one  cr  ellipsoid 
is  greater  than  0.  19874.  One  such  sphere  whos  radius  is  always  larger 
than  any  o'-  is  one  whose  radius  is: 


d.m,  =  (cr? 
rms  *  p, 


+  <-  *  "p,* 


P2 


where  (ot*  +  ap~+  dp3)  *8  t*ie  trac«  of  the  covariance  matrix  A_.  The  trace 
of  an  orthogonal  matrix  is  invariant  under  a  similarity  transformation  [lOj, 
therefore, 


drms  =  <°x,  +  <7  +  °x3  $ 

It  is  possible,  under  some  conditions,  to  get  a  conservative  estimate  of 
a  system  without  using  the  R  matrix  which  will  greatly  simplify  calculations. 
For  a  system  where  the  probability  of  finding  a  position  determination 
must  be  large  (say  0.  95  or  more)  the  drmB  statistic  can  be  a  meaningful 
tool  in  providing  a  figure-of- merit.  For  example,  assume  that  a  system 
must  provide  a  probability  of  0.97  that  a  position  determination  is 
within  10  miles  of  the  true  position  for  some  given  region  of  space.  It  is 
clear  that  a  system  whose  three  a  ellipsoid  is  contained  in  a  sphere  whose 
radius 

D(0.97)  =  3drm»  =  [(3crX|)J  +  (3ctX2)2  +  (3ctX3)2J^ 

is  less  than  10  miles,  is  a  useful  system.  Use  of  the  drmB  statistic  is 
made  whenever  system  considerations  allow  (see  reference  [ll]). 


APPLICATIONS 

1.  RHO-RHO-RHO  Satellite  System.  The  system  to  be  analyzed 
will  consist  of  two  satellites  used  as  follows:  A  ground  station  would 
transmit  a  radio  signal  to  a  satellite  which  would  relay  it  to  an  aircraft. 
The  aircraft  would  receive  and  retransmit  the  signal  and  the  aircraft 
altitude  back  to  the  original  satellite  and  also  to  a  second  satellite.  They, 
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in  turn,  would  relay  the  signals  back  to  the  ground  station.  Knowledge 
of  the  satellite  positions  and  the  total  delay  times  involved  in  sending 
and  receiving  the  signal  is  used  to  determine  the  distances  between  air¬ 
craft  and  satellites.  From  this  information  and  knowledge  of  aircraft 
altitude,  a  position  can  be  determined  [3]  and 

The  three  ranges  from  which  a  position  determination  can  be  made 
are  the  two  distances  between  the  aircraft  and  two  satellites  and  the 
distance  from  the  center  of  the  earth  to  the  aircraft,  which  must  be  found 
frcm  aircraft  altitude.  Because  the  earth  is  not  a  perfect  sphere,  an 
iterative  computational  procedure  must  be  followed  to  get  an  accurate 
absolute  fix  on  aircraft  position.  When  a  position  determination  is  to 
be  made,  a  computation  is  made  using  the  two  satellite -air craft  ranges 
and  the  aircraft  altitude  added  to  some  average  value  of  earth  radius. 

The  approximate  position  is  thus  determined  and  a  new  value  of  earth 
radius  for  the  approximate  position  is  found.  Another  computation  is 
made  using  the  two  satellite-aircraft  ranges  and  the  aircraft  altitude 
added  to  the  new  value  of  earth  radius.  An  accurate  position  determination 
is  thus  made.  If  necessary,  a  number  of  iterations  can  be  made  refining 
the  position  determination  until  an  accurate  reading  is  found.  In  practice, 
it  would  be  quite  simple  to  determine  the  accuracy  of  ^  reading  simply 
by  examining  the  result  of  two  or  three  iterations  of  the  computation. 

When  the  difference  between  results  is  negligible,  the  result  is  accurate. 

In  order  to  properly  use  the  probabilistic  system  model  developed 
in  the  previous  part  of  this  report,  it  is  necessary  to  identify  the  $,  f  , 

0  angles  of  the  T  matrix  with  a  set  of  system  angles. 


FIG.  10  A  TWO-SATELLITE  SYSTEM  SHOWING  THE  SYSTEM 
ANGLES  IDENTIFIED  AS  THE  ♦  AND  +  ANGLES  OF 
THE  T  MATRIX  OF  THE  PROBABILISTIC  SYSTEM 
MODEL 
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Fig.  10  is  a  sketch  of  a  two-satellite  RHO-RHO-RHO  system 
showing  the  system  angles  which  are  chosen  to  correspond  to  the  T 
matrix  <J>,  ^  angles  of  the  probabilistic  system  model.  They  are  the 
angles  between  the  R|,  R3  and  the  R2,  R3  axes,  respectively.  Fig.  10 
also  shows  the  subsatellite  points  Sj  and  S2  which  are  the  points  of 
intersection  of  the  earth  sphere  and  the  lines  Si  S3  and  S2  S3 ,  respectively. 
The  9  angle  of  the  T  matrix  corresponds  to  the  satellite  system  angle 
formed  by  the  tangents  of  the  great  circle  arcs  drawn  between  the  sub¬ 
satellite  point  Sj  and  PQ  and  the  subsatellite  point  S;  and  P0  as  in 
Fig.  11. 


FIG.  11  THE9  ANGLE  OF  THE  T  MATRIX  CORRESPONDS 
TO  THE  SYSTEM  ANGLE  BETWEEN  THE  TAN¬ 
GENTS  OF  THE  GREAT  CIRCLE  ARCS  DRAWN 
BETWEEN  Sj  AND  PQ  AND  BETWEEN  s'2  AND  PQ 


Fig.  11  also  shows  the  x,,  x2  axes  of  the  orthogonal  coordinates 
(X|,  x2,  X3)  into  which  errors  will  be  transformed  by  the  T  matrix.  That 
is, 


0 

esc  ^  esc  6 
0 


where  R3  and  x3  are  coincident  and  normal  to  the  earth  sphere  at  PQ  and 
where  x,  and  x2  are  orthogonal  to  each  other  and  to  x3  and  in  a  plane 
tangent  to  the  earth  at  PQ  with  x,  coincident  with  the  tangent  of  the  great 
circle  arc  5|po  at  PQ.  The  errors  in  the  (x,,  x2,  x3)  coordinate  set  are 
therefore: 


Ci,  =  Cj,  csccj)-  cjj3  cot  $ 

Cx2  =  “C*,  esc  $  cot  0  +  C*2  esc  esc  9 
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+  Cj3  cot  (J)  cot  0  -  Ci3  cot  CSC  0 


Under  the  assumption  that  for  air  traffic  control  interest  is  in  the  position 
of  an  aircraft  in  the  plane  tangent  to  the  earth  at  P0,  it  is  necessary  to 
consider  only  the  errors  in  such  tangent  plane.  For  the  two-satellite 
system,  therefore,  onlyCx  an<*  Cjj  •  the  errors  along  the  orthogonal 
axes  in  the  tangent  plane  need  to  be  Considered.  Since  the  ellipsoid  is 
a  surface,  all  plane  sections  of  which  are  ellipses,  the  loci  of  constant 
probability  error  density  in  the  plane  defined  by  the  coordinates  x,  and 
x2  are  ellipses.  They  are  called  bivariate  normal  cqui- probability  error 
density  ellipses.  The  two-dimensional  probability  density  function  for  a 
normal  distribution  in  the  (x,,  x2)  coordinate  set  is: 


-Xj)  = 

.  i 

e  7 

q£x,,Cx2) 

2 *  <y*\axi 

Xj)  = 

1  7 

r— s  .  CxrCxir 

O.-i.  1  X 

1*1  ' 

'  ^Xi 

1,  k 

e  = 

^11  ^12 
^21  ^22 

'  N 

9  jk  ~ 

1 

Tel 

axj  ff*kN  1=1  J 

1 

(j  =  1,2) 

(i,  k  =  1,2) 


If  the  errors  were  known  in  the  two-dimensional  principal  axes  coordinate 
set  (p, ,  p?)  then: 


e  ■  r  k'i) 

'jk  t  1  k  =  j  f 


(j,  k  =  1,2) 


Therefore, 
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1 

TT 


(j  =  1.2) 


where  <r 2 

Pj 


The  probability  density  function  in  the  principal  axes  set  is: 

-2 


p  £p,£p?) 


.J*,  +  CM 

-  Ss  57 


2  » 


'P, 


To  find  the  probability,  P  ,  that  a  position  determination  is  in  an  ellipse 
whose  axes  are  c  (Tp(  and  c  Opj  where  c  is  a  constant,  evaluate  the  integral 


in  polar  coordinates 


z,  =  r  cos  8 
z2  =  r  sin  © 


+  Zj  -  r? 
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let  y 


r? 

2 


dy  =  r  dr 


The  quantity  in  brackets  can  be  found  in  tables  of  Probability  and 
Statistics  books .  A  short  table  of  Pc  values  is  given  in  Table  II. 

TABLE  II 

THE  PROBABILITIES,  P^,  THAT  A  POSITION 
DETERMINATION  F\LLS  IN  A  ONE,  TWO,  TWO  AND  ONE- 
HALF  OR  THREE  SIGMA  ELLIPSE 


Table  II  shows  that  the  probability  ihat  a  position  determination  falls 
into  a  one,  two,  two  and  one-half,  and  three  sigma  ellipse  is  .  3935, 

.8647,  .9561  and  .9889,  respectively.  The  physical  axes  can  be  found 
by  using  the  R  matrix  of  the  system  model  developed  in  the  previous 
part  of  this  report.  To  simplify  computations  assume  that  a  satellite 
data  acquisition  system  should  provide  a  minimum  probability  of  .95  that 
a  position  determination  is  within  some  specified  distance  from  the  true 
position.  A  circle  which  contains  the  2.  5  sigma  ellipse  (.9561  probability) 
is  one  whose  radius  is  D  -  2,  5  drms  where' 

dr™,  =  [oj,  * 


2S 


where 


x, 


r2  - 


N 

N~  S  (C  x,)2 
i  =  i 
N 


iSr^ (C  ** 

i  =1 

Cx,=  G  J,  esc*  -  ci,  cot<J> 

«c$cot0  +  cijcsc  4>  csce 

+  ci cot  $ cot  b  *  ci  cot  ^ csc  ® 

"3  *3 

N 

°x  =  (Gi  esc*  -  ei  cot*)2 

*1  N  ^  "1  *3 

a2  =  crj^  esc-  *  +  cr*3  cot2* 

sinceG^i  G  J3  are  independent  random  variables,  i.e., 

N 

"  ^-Ci,  Gi3  CSC  *  cot  *  =  0 


1  =1 


Similarly, 


°X2=  Oj,2  CSC2*  cot20  +  CSC2  V'  CSC2  0  +  (Tj3  (cot  *  cot  0  -  cot  ^  CSC  0  )2 


D  =  2.  5  d 
D 


rms 


2.5  csc2*  +  <r2^cot2$ 

+  <7^  CSC2  *C0t2  0  +  (Tj^CSC*^  CSC20 

+  (T,2  cot2  *  cot2  0  +  oj  cot2  ^  CSC2  0 

3  3 


-  2  (T 


23  cot  *  cot  0  cot  +  csc0j7 


D 


=  2.5  £<y csc2*  csc20 

+  <t|2  CSC2  ^  CSC2  0 


-L  cr: 


(cot2*  csc2  0+  cot2  ^  csc20  -  2  cot*  cot  ^  cot  0csc0)J 
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FIG.  12  PQ  IS  THE  TRUE  POSITION,  Pi  THE  MEASURED 
POSITION  DUE  TO  ERROR  IN  MEASURING  THE 
STATION  POSITION  ALONE 


From  Fig.  12: 

=  |C!„,|  co.  1-  *  y 

y  =  (Ri  +  j cj0| j co»  > )  ( 1  -  co»  e) 

Therefore, 

Cj,  =  |cj0(j  cos  y  +  (R,  +  |  ^so,|  cos  >)  (1  -  cos©) 
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whereGsoiis  the  i'th  error  vector  to  be  determined  for  the  particular 
system  used  to  measure  the  position  of  station  one.  For  small  0,  Csi 
j  C  soil  cos  y  ,  which  is  the  component  of  the  error  vector  on  R,.  For  the 
satellite  system  under  consideration,  where  the  ranges  are  thousands  of 
miles  and  the  errors  are  of  the  order  of  a  few  miles  (a  conservative 
estimate),  0  is  exceedingly  small  and€|t  ,  €  ,  £  j3  can  be  taken  as 

the  components  of  the  error  vectors  ^j0|.  ^soj*  ^$03  along  R|»  Rj.  and 
R3,  respectively.  The  errors  in  the  i'th  measurements  on  the  ranges 


are: 

3  d,+  ci, 

=  d,+  d. 

^3T 

•  dj+  d, 

Therefore: 

=  ®2,  +  o-l, 

=  *  al , 

-  « 2 ,  -  o-l, 

Since  the  correlations 

^R.  S- 
J  J 

1 

N 

_ y  e 

N  <7t 

"J 

,  <rs  2-, 

1  5I  1  =  1 

Therefore: 

1 


(j  =  1,2,3) 


=  2.5  £  (t»)t  csc2<t>csc20 

+  a^2T  esc2  esc2  0 

+  (cot2$  csc20  +  cot2  ^  esc2  0  —  2  cot$  cot  ^  cot  0csc0)j^ 


In  order  to  get  a  mapping  of  error  isograms  (constant  D  contours)  on  the 
earth  it  is  necessary  to  find  functional  relationships  between  the  $,  ^  , 
and  0  parameters  to  some  set  of  earth- related  parameters. 
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FIG.  13  THE  SYSTEM  ANGLE  <& IS  RELATED  TO  THE 
EARTH  ANGLE  0 

From  Fig.  13: 

a3  =  Rj  +  Rj  -  2  R,  R3  cos  <X> 

R2  =  a3  L  Rj  -  2  a  R3cos  0 
Eliminating  RI  and  solving  for  <J>  gives: 


cos  = 


a  cos  0  -  R3 
(a2  +  R3  -  2aRjCos  0$ 


FIG.  14  THE  SYSTEM  ANGLE^IS  RELATED  TO  THE 
EARTH  ANGLE  > 


From  Fig.  14: 

b2  =  R22  +  Rl  -  2  R2  R3  cos  * 
Rj  =  b2  +  R|  -  2  b  Rj  cos  J 


(6) 
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Eliminating  R2  and  solving  for  p  gives: 


b  cos  7  -  R3 

cob  i  -  (b2  +  R2  _  2  b  Rj  CQS  y ) 

With  the  help  of  Fig.  15  it  is  possible  to  transform  from  0,  y  ,  0  to 
L,  x,,  X2  a  set  of  earth  angles  which  become  Latitude  and  Longitude  angles 
when  the  two  satellite  stations  are  directly  over  the  equator. 


FIG.  15  THE  ANGLES  0,  7  RELATED  TO  THE 
L,  X,,  \2  COORDINATE  SET 

From  Fig.  15  using  spherical  trigonometry: 

ccs  (X)  +  Xj)  =  cos  0  cos  y  *  sin  0  sin  7  cos  0 

cos  0  -  cos  L  cos  X i 

cos  7  =  cos  L  cos  X2 

Using  equations  (6)  and  (7)  and  the  equations  derived  from  Fig.  15: 
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cos  0 


cos  (XjJ-_Xj)  —  cos  L  cosX)  cos  L  cos  X2 
pi  -  cos2  L  cos 2  X,)  ( 1  -  cos2  L  cosJX2)Jl5 

a  cos  L  cos  X,  —  R3 _ 

(a2  R^  -  2  a  R3  cos  L  cos  X,)^ 

cos  ^  _  b  cos  L  cos  X  j  -  R3 _ 

(b2  +  R|  -  2  b  R3  cos  L  cosXTjH” 

Summarizing: 

D  =  2.  5  ^ o^T  esc2  4>  esc2  0 

+  <rJ2T  csc2  ^  cs<^0 

+  cr2^T  (cot2<J>csc20  +  cot2^csc^0  -2  cot<£> cot^ cot0csc0)J$ 

where: 


cos  0  = 


cos  (X,  +  X  j  -  cos2  L  cos  X2  cos  X) _ 

(1  -  cos 2  L  cos2  X,)  ( 1  -  cos2  IT  cos7X2V'2 


cos<D=  g--  c£-s  -  L  -co-8-*  1  --  R? _ _ _ 

(a7  +  R  l  -  2  a  R3  cos  L>  cos  X,)*} 

.  _  _  b  cos  L  cos  X2  -  R3 _ 

COS  (b2  +  R^  ”  2  b  Rj  cos  L  cos  Xj)  ^ 

where  a  and  b  are  the  distances  between  the  earth  center  and  the  satellites 
S,  and  S21  respectively.  A  computer  program  was  written  for  a  two- 
satellite  system  where  the  satellites  were  spaced  at  angles  ranging  from 
40  to  70  degrees  on  the  equator  in  synchronous  altitude  orbits.  To  map 
constant  D  values  (error  isograms),  it  is  necessary  to  make  an  estimate 
on  the  and  Os ,  >  O'sj*  ffS3‘  References  [3]  ,  (4  ]  and  [5]  include 

detailed  analyses  and  studies  from  which  the  values  can  be  estimated. 

It  is  assumed  that  the  satellite  ranges  can  be  measured  to  within  300  feet 
(one -sigma  value)  and  the  aircraft  altitude  can  likewise  be  measured  to 
within  300  feet  (one -sigma  value).  Since  more  sophisticated  equipment 
can  be  used  to  measure  the  satellite  positions,  it  is  assumed  that  0‘s,, 

CTjj  are  100  feet  (one-sigma  value).  0*3,  the  measure  of  the  position  of 
the  earth  center  is  assumed  known  to  such  a  close  value  that  it  is  taken 
as  zero.  That  is,  it  is  assumed  that: 


^3 


0*3  =  300  feet 
100  feet 


OSj 

0 
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Using  these  values,  mappings  of  constant  D  values  in  nautical  miles  are 
shown  for  the  North  Atlantic  Ocean  Area  in  Figs.  17  to  20  for  angular 
satellite  spacings  of  40,  50,  60,  and  70  degrees.  Because  the  values  of 
assumed  standard  deviations  are  rather  optimistic  and  are  based  upon 
an  attempt  to  minimize  them,  the  set  of  mappings  in  Figs.  21  to  24  are 
based  upon  a  more  conservative  set  of  assumed  standard  deviations  as 
follows: 


<r.  =  cr,  =  <r»  =•  900  feet 

*1  *2  *3 

cr j  =  <Tj2  =  300  feet 

Or  =  0 

*3 

When  a  position  determination  is  made,  the  probability  is  greater  than 
.95  that  it  is  within  D  nautical  miles  of  the  true  position. 

Figs.  17  to  24  show  that  excellent  results  are  attainable  for  the 
two-satellite  system  over  a  large  portion  of  the  North  Atlantic  Ocean 
including  the  heavily  travelled  air  route  between  Newfoundland  and  Ireland. 
Far  north  regions  are  excluded  and  large  D  values  are  indicated  near  the 
equator.  The  equation  for  D,  in  fact,  indicates  that  at  a  position  on  the 
equator,  D  is  infinite.  The  reading  is  physically  absurd  and  is  due  to  the 
first  order  approximation  used  in  constructing  the  probabilistic  model  in 
the  General  Analysis.  In  finding  the  error  vector  for  a  particular  mea¬ 
surement  it  was  assumed  that  for  large  spheres,  small  portions  of 
spherical  surfaces  could  be  considered  planes.  The  assumption  leads  to 
the  singularity  in  D  at  the  equator.  To  get  an  upper  bound  on  the  error 
that  can  be  expected  at  the  equator,  it  is  necessary  to  apply  a  second 
order  approximation.  It  is  apparent  from  geometrical  considerations 
that  the  subsatellite  points  s!  and  Sj  shown  in  Fig.  16  will  be  measured 


FIG.  16  THE  SUBSATELLITE  POINTS  S;,  SJ  ON  THE 
EQUATOR  ARE  MEASURED  WITH  THE  LEAST 
PRECISION 
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The  probability  is  greater  than  .95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  cr^  =  <7*5  =  a*  =  300  feet,  0$  = 

Oj  =  100  feet,  Os3  =  0  feet. 


FIG  18  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  0°  WEST  LONGITUDE 
AND  60°  WEST  LONGITUDE  -  I 
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The  probability  is  greater  than  .95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  *ext  for  definitions):  o„)  =  cr^  =  ct*3  =  300  feet,  crS|  = 
<TSj  =  100  feet,  c-S;j  =  0  feet. 


FIG.  19  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  5°  WEST  LONGITUDE 
AND  55°  WEST  LONGITUDE  -  I 


The  probability  is  greater  than  .95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  cm  =  a;  =  <r,  =  300  feet,  or,  = 

*i  *2  3  *1 

<TSj  =  100  feet,  a$3  =  0  feet. 


FIG.  20  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  I’HE  NORTH 
ATLANTIC  AREA  MUTUALLY  VISIBLE  TO  SATELLITES 
STATIONED  ABOVE  10°  WEST  LONGITUDE  AND  50°  WEST 
LONGITUDE  -  I 
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The  probability  is  greater  than  .  95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  <T,|  -  <Tnj  =  =  900  feet, 

o^-  300  feet,  < rS;)  -  0  feet. 


FIG.  21  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  5°  EAST  LONGITUDE 
AND  65°  WEST  LONGITUDE  -  II 


The  probability  ie  greater  than  95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  <7^  =  °*  s  -  °i»3  "  900  feet,  <rS)  - 
Oj^  -  300  feet,  oSj  =  0  feet. 


FIG.  22  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  0°  WEST  LONGITUDE 
AND  60°  WEST  LONGITUDE  -  II 
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The  probability  is  greater  than  95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  O’.  =  crB  =  o.  =  900  feet.  <x  = 

si  si  J  | 

<t.  -  300  feet,  a,  =  0  feet. 

*2  >3 


FIG.  23  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  5°  WEST  LONGITUDE 
AND  55°  WEST  LONGITUDE  -  II 


The  probability  is  greater  than  .95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  input  standard  deviations 
as  follows  (see  text  for  definitions):  OJ  =  =  <T„  =  900  feet,  <rs  = 

<rSj  =  300  feet,  c r$j  =  0  feet.  i  a  ) 


FIG.  24  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  MUTUALLY  VISIBLE  TO 
SATELLITES  STATIONED  ABOVE  10°  WEST  LONGITUDE 
AND  50°  WEST  LONGITUDE  -  II 
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Since  the  earth's  sphere  is  much  smaller  than  the  spheres  generated 
by  either  R,  or  R>,  a  useful  approximation  can  be  made  by  assuming  that 
the  portion  of  the  sphere  R|  =  constant,  near  S'|t  is  a  plane.  Fig.  25  shows 
the  approximation. 


FIG.  25  THE  REGION  ABOUT  THE  POINT  S',  SHOWING 
MEASUREMENT  DEVIATIONS  A  R,  AND  A  R3 

In  Fig.  25  A  R,,  A  R3  are  deviations  in  the  measures  of  R,  and  R3. 
respectively. 

For  <T*(  =  300  feet  and  oS)  =  100  feet 
A«,  =  2.2  Vcr*,  +  ai; 

A*,  =  2.2  V  ( 3 00)2  +  (Too?1 

As,  =  695  feet 

The  probability  is  .972  that  the  measurement  of  R,  is  within  695  feet  of 

the  true  value  of  R,. 

For  <Ts3  =  300  feet  and  crS3  =  0 

Ar3  =  2.2  ^ <7s3  +  o’*, 

As3  =  660  feet 

The  probability  is  .972  that  the  measurement  on  Rj  is  within  660  feet  of 

the  true  value  of  Rj.  The  probability  is  .95  that  the  measurement  of 
position  at  Sj  is  within  A  Rt  and  within  A  Rj  of  the  true  values  of  R,  and 
R3.  The  probability  is  greater  than  .95  that  the  measurement  is  within 
the  distance  d  of  the  true  position  where  d  is  shown  in  Fig.  26. 
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FIG.  26  THE  DISTANCE  d  SHOWS  THE  LARGEST 
LINEAR  DISTANCE  PERPENDICULAR  TO 
R3  (AND  R,) 

From  Fig.  26: 


cos  0  = 


r3  -ar, 

R3  +  AR| 


sin  0  = 


d  = 
d  = 
d  = 
d  = 
d  = 
d  = 
d  = 


d 


R3  +  AR, 

(R3  +  AR,)  (1  -  co8J0),y2 


(R,  +ARI)[1-  (fj- 

[(R3  +  A  R,)J-  (R3  - 
[r?  +  2  R3  A  R,  +  A  R?  - 


-  AR, 

+  AR, 

AR,)’! 

R3+  2  R3AR, 


2  RsAR,1 
2  V*000  X  .  132' 


39.8  nautical  miles 


For  7»,=  300  feet  and  7$,  =  100  feet  the  probability  is  greater  than  .95 
that  the  measured  position  is  within  39.8  nautical  miles  of  the  true  posi¬ 
tion  when  the  true  position  is  the  subsatellite  point  S|.  It  should  be 
emphasized  that  the  value  is  an  upper  bound  and  that  real  values  can  be 
determined  only  by  experiment. 
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2.  RHO- THETA-THETA  Satellite  System.  The  system  consists 
of  a  single  satellite  with  orthogonal  interferometer  arms  and  a  distance - 
measuring  capability.  From  a  knowledge  of  the  angles  between  the  air¬ 
craft  and  each  interferometer  arm  and  the  distance  between  satellite  and 
aircraft,  a  position  can  be  determined  (for  details  see  reference  [f]). 


z 


y 


FIG.  27  AN  INTERFEROMETER  SYSTEM  MEASURING 
THE  ANGLES  A  AND  B  AND  THE  DISTANCE  ? 

Fig.  27  is  a  sketch  of  an  orthogonal  interferometer  system  with 
the  interferometer  arms  along  the  x  and  y  axes  and  the  z  axis  orthogonal 
to  both  arms.  The  errors  shown  are  linear  errors  with  € $  (A  A)i 
and  =  9(AB)*  where  (A  A)*  and  (AB)*  are  angular  errors  in  A  and  B, 
respectively.  The  coordinate  set  of  the  error  components  is  not  orthogonal 
making  it  necessary  to  apply  the  T  matrix  of  the  probabilistic  system 
model  developed  in  the  General  Theory.  To  use  the  T  matrix,  it  is  nec- 
cessary  to  identify  system  angles  with  angles  of  the  T  matrix  as  follows: 
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a.  Identify  <t>  as  the  angle  between  e^and  =  90° 

b.  Identify^  as  the  angle  betweens^and  =  90° 

c.  Since  <t>  and  4/  are  both  90°,  the  angle  0  is  the  angle 
betweenGiandej.  The  T  matrix, 


T  = 


CSC  $ 

-  CSC  <J>  cot  0 
0 


0 

CSC  4*  esc  0 
0 


-  cot  <t> 

cot  <X>  cot0  -  cot  4*  esc 0 
1 


is  therefore: 


T  = 


1 

-cot  0 
0 


0  0 

esc  0  0 

0  1 


thus: 


and 


1  0  o" 

ef 

-  cot  0  CSC  0  0 

e  A 

°  0  1  _ 

St 

= 

=  -cj  cot  0  +  f"  i  CSC  0 


are  the  error  components  in  the  orthogonal  coordinate  set  (x,,  Xj,  x3) 
where  xi  is  along  B,  x3  is  along  ?and  x3  is  orthogonal  to  both  X|  and  x3 
as  shown  in  Fig.  28. 


44 


FIG.  28  THE  ORTHOGONAL  SET  (x„  x}(  x,)  RELATED 
TO  THE  NON -ORTHOGONAL  SET  (?,  A,  B) 

In  Fig.  27  the  earth's  center  is  located  somewhere  along  the  z  axis, 
and  the  earth's  sphere  intersects  the  spherical  segment  of  Fig.  27  in  the 
arc  EE'.  It  is  evident  that  none  of  the  orthogonal  error  components 

is  in  the  tangent  plane  of  the  earth's  surface.  To  get  cjt,  = 

G  £  in  the  tangent  plane  of  the  earth's  surface,  the  system  must  be  rotated 
counterclockwise  about  the  £  axis  through  the  angle  a  shown  in  Fig.  27. 

Then: 


COS 

a 

sin  a 

0  ' 

1 

n 

X  — 

i 

CV, 

= 

-sin 

a 

cos  a 

0 

i 

€x2 

-CV, 

0 

0 

1 

ex 

L  X3  J 

i 

CYi 

= 

< 

cos 

a  + 

C  h.  sin 
xi 

a 

CY, 

= 

Gx, 

s  in 

a  + 

£  i  cos 

Xj 

a 

eVj 

= 

X3 

where  C.  ^  (  is  in  the  tangent  plane  of  the  earth  as  shown  in  Fig.  29. 


WITH  G  1  IN  THE  TANGEN  T  PLANE  OF  THE 
EARTH  ' 

From  Fig.  29,  c  y2  can  be  placed  into  the  tangent  plane  of  the 
earth  by  rotating  clockwise  about  the  y,  axis  through  the  angle  0  .  Thus: 


and 

c  y. 

ey}  cos  0  -  sin  0 

G  X  sin  0  +  C  ^  cos  0 

where 

and  are  in  the  tangent  plane  of  the  earth. 


46 


=  c 

1  cos 

y? 

0 

- 

cy> 

sin  0 

< 

=  c 

i 

cos 

X| 

a 

+ 

Cx 

x3 

sin  a 

=  ( 

-  c*, 

sin 

a 

+  c 

t  cos 
x2 

1 


C'  sin  0 


=  Cl  (cos  a  -  cot  0  sin  a)  +  Ca  C8C  &  sin  a 
=  Cl  (  -  sin  a  cos  0  -  cot  0  cos  a  cos  0  ) 

+  Ca  esc  0  cos  a  cos  0 
-  C^>  sin  0 


Similarly, 


1  =1 


=  a\  (cos  a  -  cot  0  sin  a  )2  +  ®2  (esc  0  sin  a  )* 


(3  = 


(J  j  ( —  sin  a  -  cos  a  cot  0  )2  cos2  0 
a\  (esc  0  cos  a  cos  0  )2 


+  sin2  0 


rms 


+  o 


d2 

rms 


a 2  £(cos  a  -  sin  a  cot  0  f  +  (-  sin  a  -  cos  a 

cot  0)s  cos2  dl 

+  <J2  |"csc20  sin2a+  csc20  cos2<*  cos2dj 

+  a £  |s  in  2 

For  a  measurement  to  have  a  probability  greater  than  .95  of  being 
within  some  distance  D  of  the  true  position, 


2. 5  d 


rms 
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Then: 

D  =  2.  5  |  <x2  £(cos  a  -  sin  a  cot  0  )*  +  (  -  sin  a  -  cos  a  cot  0  )J  cos2/? 
+  <5  2  £csc2  0  sin2a  +  csc20  cos2a  cos2  $  J 
+  a|£sin20jj&  (8) 

To  get  a  mapping  of  error  isograms  (constant  D  curves)  on  the 
earth's  surface,  it  is  necessary  to  relate  the  angles  0,  a>  fi  to  a  system 
of  earth  angles.  From  Fig.  27,  using  spherical  trigonometry: 

cos  90°  =  cos  A  cos  B  +  sin  A  sin  B  cos  6 

cos  O  =  -  cos  A  cos  B 
sin  A  sin  B 


cos  A  =  cos  P  co<®  C 


sin  A 


=  V~1  -  cos2  C  eoe2  P1 


cos  B  a  cos  (90°  -  P)  cos  C 

cos  B  =  sin  P  cos  C 

sin  B  =  -  sin2  P  cos2 


Therefore: 

cos2  C  sin  P  cos  P 

cos  0  =  -  - -  -  - —  - - -  - - - - -  -  (9) 

^  1  -  cos2  C  cos2  P1  ^  1  -  cos2  C  sin2  P 


From  Fig.  27  the  maximum  value  of  C  is  90°.  The  minimum  value  of  C 
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From  Fig.  30  with  /i  -  90° 


a2 

-  R2  +  K>7 

P2 

=  R2  +  a2 

-  2  a 

R 

cos  y 

a2  -  R2 

=  R2  +  a2 

-  2  a 

R 

cos  y 

R 

cos  y  =  ~~ 

cl 

For  a  satellite  in  a  synchronous  orbit: 

4000 

cos  7  *  26300 
cos  y  «  .151 

y  «  81° 

From  Fig.  30: 

90°  -  Cmin  =  180  -  90  -  y 

cm,„  =  y 

Cmin  «  81° 

COS  Cm,„  w  .  151 
COS2  Cmi n  «  .023 

To  a  very  good  approximation  let  cos2  C  =  0  in  the  denominator  of 
equation  (9).  That  is,  in 

_  cos2  C  sin  P  cos  P 

cos  W  -  -  -  -  . -  1  ■  1  ■  ‘  - ; 

•^1  -  cos2  C  sin2  P  yl  —  cos2  C  cos2  P 
let  cos  C2  =  0  in  the  denominator.  Then: 

cos 6  =  -  .023  sin  P  cos  P 

The  maximum  value  of  cos  0  is  given  for  P  =  45°.  Therefore, 
cos  0mo„  =  -  .  023  (  ) 
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cos  H 


.  0115 


"'ox 


90°  40’ 


Therefore: 


90°  >  0  >  90°  40' 

That  is,  the  value  of0is  very  close  to  90°  over  the  entire  surface  of  the 
earth  when  the  satellite  is  in  a  synchronous  orbit. 

In  equation  (8)  for  D,  the  multipliers  of  cotOand  esc  0  are  sin  a 
and  cos  a .  Since  the  maximum  values  of  these  functions  is  unity  to  a  very 
good  approximation,  let 

cot  0=0 

esc  0=1 


Then: 


D  =2.5  |  crB2  £  cos2  a  +  sin2  a  cos2  /(jj 
sin2a  +  cos2  a  cos* 
sin7 


+  °A2 


+  cri 


Because  the  angles  A  and  B  are  both  large  over  the  entire  earth  surface, 
assume  that  Oa  =  0*  =  O  then, 


D 


2.5 


1  +  cos2  0 


]+ 


(10) 


To  a  very  good  approximation  the  D  values  are  independent  of  the 
angle  P  and  the  loci  of  constant  D  «.re  concentric  circles  on  the  earth's 
sphere  centered  at  the  subsatellite  point  on  the  earth's  sphere. 


Reference  [5]  contains  an  extensive  analysis  of  factors  which  affect 
<T  and  Ofa  in  order  to  arrive  at  some  reasonable  estimate  for  those  values. 
The  values  assumed  in  reference  [5]  are  approximately: 


a  =  £  x  .  03  milliradians 


=  300  feet 

Since  the  minimum  value  of  £  is  22,  300  miles, 
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<x  =  2?,  300  X  5,280  X  .03  X  10“  3  feet 
a  «  3,  532  feet  (a  minimum  value) 

The  maximum  influence  of  a ^  on  D  in  equation  (10)  is  exerted  at  £  =  90°. 
At  that  point: 

D  =  2.  5  (<y2  +  cr2  yi 

Since  cr^o^  the  influence  of  c^>  on  the  system  is  very  slight  and  system 
precision  is  basically  a  function  of  the  angular  errors  of  the  interferometer 
arms.  Assuming  is  negligible, 

D  =  2.  5  cr  (1  +  cos2  0)%  (11) 

where 

<r  =  .03  x  io~ 3  e 

_ „  a  R  -  a  cos  7 

cos  p  =  - — — - 

e 

£  =  (a2  +  R2  -  2  a  R  cos  7  )* 

R,  a,  and  7  are  shown  in  Fig.  30.  Constant  D  values  are  plotted  in  Fig. 

31.  The  probability  is  greater  than  .95  that  the  measured  position  is 
within  D  nautical  miles  (to  a  very  good  approximation)  of  the  true  position 
where  the  satellite  is  in  a  synchronous  orbit  at  30°  West  Longitude. 

There  exists  no  known  experimental  data  on  the  performance 
of  an  interferometer  in  space,  such  as  the  system  proposed  here, 
or  even  of  a  spaceborne  platform  which  measures  a  single  angle.  How¬ 
ever,  some  experimental  data  exists  on  the  accuracy  of  a  ground-based 
interferometer  system  in  use  at  NAFEC  (National  Aviation  Facilities 
Experimental  Center)  of  the  Federal  Aviation  Agency  located  at  Atlantic 
City,  New  Jersey.  MOPTAR  (Multi-Object  Phase  Tracking  and  Ranging) 
built  by  the  Cubic  Corporation  uses  interferometer  arms  which  are  125 
feet  long.  The  angular  accuracy  claimed  and  corroborated  by  testing  at 
NAFtuC  for  absolute  position  determination  is  200  parts  per  million  in 
direction  cosine  (reference  [12J).*  For  angles  near  90°,  200  parts  per 
million  in  direction  cosine  is  equivalent  to  0.  2  milliradians .  The  estimate 
of  .03  n-illiradians  for  a  200-foot  long  spaceborne  interferometer  seems 
very  optimistic.  Assuming  that  the  tolerances  could  be  held  on  a  space¬ 
borne  125-foot  interferometer  system  and  the  angular  deviation  for  the 
system  is  0.  2  milliradians,  the  equation  (11)  gives  constant  D  values  as 
shown  in  Fig.  32. 

*  An  integral  part  of  the  measuring  technique  consisted  of  calibration  of 
the  system  against  fixed,  surveyed  targets  before  and  after  each  test  run. 
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The  probability  is  greater  than  .95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  an  input  angular  standard 
deviation  of  .  03  millira  ms  . 


FIG.  31  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  VISIBLE  TO  A  SATELLITE 
STATIONED  ABOVE  30°  WEST  LONGITUDE 
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The  probability  is  greater  than  .  95  that  a  measured  position  is  within 
D  nautical  miles  of  the  true  position  assuming  an  input  angular  standard 
deviation  of  .  2  milliradians  (experimental  value  of  the  NAFEC  MOPTAR 
ins  tallation). 


FIG.  32  ERROR  ISOGRAMS  (CONSTANT  D  VALUES)  FOR  THE 
NORTH  ATLANTIC  AREA  VISIBLE  TO  A  SATELLITE 
STATIONED  ABOVE  30°  WEST  LONGITUDE  -  125-FOO,p 
INTERFEROMETER 
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DISCUSSION 


The  position  determination  probabilistic  system  model  developed 
in  this  report  is  a  useful  computational  tool  with  the  proper  identification 
of  the  model  angles  ($,  ,  9  angles  of  the  T  matrix)  to  the  angles  of  any 

particular  position  determination  system.  Once  the  identification  is 
made,  the  use  of  the  model  leads  to  a  solution  in  a  straightforward  manner 
and  results  are  presented  in  the  form  of  the  physical  size  and  orientation 
of  the  trivariate  normal  equi-probability  error  density  ellipsoid.  The 
model  was  but  partially  used  in  the  Applications  section  of  this  report, 
because  the  nature  of  the  solution  sought  allowed  the  use  of  the  drms  error 
statistic  which  eased  the  computational  problem.  It  should  be  empha¬ 
sized,  however,  that  the  model  docs  include  the  mathematical  method¬ 
ology  necessary  for  getting  the  more  exact  physical  picture  of  position 
determination  system  errors  in  the  form  of  the  trivariate  normal  equi- 
probability  error  density  ellipsoid  and  its  orientation  in  three-dimensional 
space,  when  system  considerations  warrant  the  computational  work  in¬ 
volved.  The  model  makes  use  of  a  first-order  approximation  in  the 
system  geometry  by  assuming  that  for  small  distances  from  the  true 
position  quadric  surfaces  of  a  position  determination  system  are  planes. 

In  applying  the  model  to  the  two-satellite  RHO-RHO-RHO  system  a 
singularity  was  found  at  the  equator  (i.e.  ,  D  =  infinity).  A  second-order 
approximation  was  required  to  get  an  idea  of  system  quality  at  the 
singularity.  The  mapping  of  constant  D  values  for  areas  other  than  near 
the  equator  shows  that  the  mot  ^1  gives  a  good  approximation  to  the  system 
errors.  Because  the  one-satellite  RHO- THETA- THETA  system  quadric 
surfaces  are  cones  with  very  large  apex  angles  (162°  to  180°),  the  diverg¬ 
ence  of  a  plane  surface  from  the  true  conical  surface  is  extremely  slight. 

A  simple  calculation  shows  that  for  D  values  less  than  25  miles,  the  plane 
does  not  diverge  from  the  cone  by  more  than  25  feet.  The  application  of 
the  model  to  the  one-satellite  system  gives  an  excellent  approximation  to 
system  errors  with  no  singularities. 

A  report  of  the  Requirements  Panel  of  the  Joint  Navigation  Sate^ite 
Committee  dated  February  1965,  recommends  aircraft  position  determin¬ 
ation  for  a  controller  to  an  accuracy  of  5  nautical  miles  in  1965,  improving 
to  3  nautical  miles  in  1975  and  to  1  nautical  mile  by  2000. 

The  predicted  D  values  for  the  RHO-RHO-RHO  technique  indicate 
that  operationally  practical  values  of  positional  accuracy  can  be  achieved 
for  attainable  instrumental  errors  in  the  system.  In  the  heavy  traffic 
density  area  between  Newfoundland  and  Ireland,  excellent  results  can  be 
expected  even  for  the  more  conservative  estimates  of  input  standard 
deviations.  Near  the  equator,  where  the  confidence  level  of  the  system  is 
low.  the  traffic  density  is  also  low  and  operational  solutions  to  the  air 
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traffic  problem  will  suffice  for  a  longer  time  than  the  operational  life¬ 
time  of  a  first-generation  satellite  system.  The  practical  aspects  of  the 
system  are  such  that  little  research  and  development  is  required  for  any 
needed  technical  innovations  or  problems.  Indeed,  the  most  difficult 
problems  for  system  implementation  for  the  two-satellite  system  would 
be  economic  and  administrative  rathen  than  technical.  As  with  any  new 
system,  some  learning  and  experience  would  be  necessary  to  solve  oper¬ 
ational  problems,  but  the  obstacles  seem  to  be  very  small  and  should  not 
limit  the  efficient  use  of  the  system.  Because  of  the  excellent  results  which 
can  be  expected  in  the  heavy  traffic  density  area,  and  the  system  simplicity, 
as  well  as  the  simplicity  of  the  satellites,  a  two-satellite  system  would  be 
a  good  choice  for  the  first-generation  communication- surveillance  satel¬ 
lite  system  for  the  North  Atlantic  Air  Traffic  Control  Area.  The  choice  of 
satellite  separation  would  depend  upon  functional  interface  between  the 
communication  subsystem  and  the  surveillance  subsystem  and  upon  the 
emphasis  placed  upon  the  possibility  of  system  expansion  for  worldwide 
coverage,  since  the  variation  in  system  precision  is  not  very  dependent 
upon  the  satellite  separation.  Future  studies  of  the  Systems  Research 
and  Development  Service  on  satellite  systems  should  be  largely 
concerned  with  the  solution  of  the  problem  of  satellite  separation  as  well  as 
other  facets  of  the  communication-surveillance  satellite  system  design. 

If  the  input  angular  standard  deviation  of  .  03  milliradians  can  be 
achieved  for  the  one-satellite  RHO- THETA-THETA  technique,  the  analysis 
indicates  that  operationally  practical  values  of  positional  accuracy  should 
be  achieved.  Available  experimental  evidence  on  the  performance  of  a 
ground-based  interferometer  (the  MOPTAR  installation  at  NAFEC)  indicates 
that  the  estimate  is  optimistic.  A  125-foot  interferometer  satellite  system 
with  the  precision  of  the  ground-based  interferometer  at  NAFEC  gives  the 
somewhat  marginal  D  values  of  Fig.  32.  Whether  or  not  a  significant  in¬ 
crease  in  system  precision  can  be  attained  by  longer  interferometer  arms 
in  the  hostile  space  environment  is  debatable.  Estimates  could  be  in  error 
by  an  order  of  magnitude  or  more.  The  satellite  system  ir  simple,  con¬ 
sisting  of  a  single  satellite  between  30°  and  40°  West  Longitude  but  the 
satellite  itself  requires  considerable  research  and  development.  The  two 
most  serious  problems  are  satellite  attitude  stabilization  and  automatic 
boom  construction  with  antenna  sensors  stragetically  distributed  along  the 
interferometer  arm  to  obviate  ambiguities.  The  problems  are  not  isolated 
from  one  another.  Attempting  to  increase  boom  length  to  increase  system 
precision  exagerates  the  attitude  stabilization  problem.  The  National 
Aeronautics  and  Space  Administration  is  planning  three  experiments  on 
gravity  gradient  stabilization  with  two  of  them  at  synchronous  orbit  altitude. 
The  ATS  (Advanced  Technological  Satellites)  series  is  planned  to  run  .'rom 
late  1966  to  late  1968.  The  research  and  development  required  on  the  tech¬ 
niques  used  by  an  interferometer -type  satellite  system,  indicates  that  from 
the  point  of  view  of  operational  traffic  control  satellites,  the  system  should 
be  considered  a  second -gene ration  system  if  a  reasonable  system  precision 
can  be  achieved  in  the  future. 
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CONCLUSIONS 


It  is  concluded  that: 

1.  A  determination  of  theoretical  system  quality  of  position 
determination  systems  can  be  made  by  the  application  of  the  mathematical 
model  constructed  in  this  report. 

2.  Error  isograms  for  the  two-satellite  surveillance  systems 
predict  excellent  results  for  determining  aircraft  position  in  the  heavy 
traffic  area  over  the  North  Atlantic  and  less  precise  measurements 

in  a  narrow  band  near  the  equator. 

3.  Experimental  evidence  on  the  NAFEC  MOPTAR  (an  interferometer 
with  orthogonal  arms)  does  not  support  the  most  optimistic  error  estimates 
for  the  interferometer  technique.  The  research  and  development  still 
required  on  the  techniques  used  by  an  interferometer -type  satellite  system, 
indicate  that  attainment  of  the  requisite  precision  by  use  of  the  one- 
satellite  technique  will  require  more  time  and  effort  that  by  use  of  the 

two  satellite  multi-rho  technique. 

RECOMMENDATIONS 


It  is  recommended  that: 

1.  The  probabilistic  position  determination  system  model 
constructed  in  this  report  be  further  exploited  by  the  Systems  Research 
and  Development  Service  for  system  studies  as  required  in  the  area  of 
air  traffic  surveillance  and  navigation. 

2.  The  Systems  Research  and  Development  Service  explore  in 
more  detail  the  two-satellite  multi-rho  surveillance  technique  to  (a) 
determine  the  feasibility  of  adapting  the  technique  for  u^e  with  a  VHF 
communication  satellite  and  (b)  formulate  a  system,  utilizing  the  technique, 
which  can  be  subjected  to  rigorous  technical,  operational  and  economic 
analyses  to  determine  its  utility  in  the  North  Atlantic  Air  Traffic  Control 
operation. 

3.  The  Systems  Research  and  Development  Service  continue 
in-house  study  of  the  interferometer -type  satellite  technique  and  also 
provide  inputs  to  experimentation  by  NASA  and  others.  Consideration 
of  the  interferometer  technique  and  comparison  with  other  techniques 
should  be  continued,  since  the  technique,  if  proven,  could  offer  certain 
advantages  in  a  later-generation  satellite  surveillance  system. 
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